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Abstract 

We examine a new optimization problem formulated in the tropical 
mathematics setting as a further extension of certain known problems. 

The problem is to minimize a nonlinear objective function, which is de¬ 
fined on vectors over an idempotent semifield by using multiplicative 
conjugate transposition, subject to inequality constraints. As com¬ 
pared to the known problems, the new one has a more general objective 
function and additional constraints. We provide a complete solution in 
an explicit form to the problem by using an approach that introduces 
an auxiliary variable to represent the values of the objective function, 
and then reduces the initial problem to a parametrized vector inequal¬ 
ity. The minimum of the objective function is evaluated by applying 
the existence conditions for the solution of this inequality. A complete 
solution to the problem is given by solving the parametrized inequality, 
provided the parameter is set to the minimum value. As a consequence, 
we obtain solutions to new special cases of the general problem. To 
illustrate the application of the results, we solve a real-world problem 
drawn from time-constrained project scheduling, and offer a represen¬ 
tative numerical example. 

Key- Words: tropical mathematics, idempotent semifield, con¬ 
strained optimization, complete solution, time-constrained project schedul¬ 
ing. 
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1 Introduction 

Tropical optimization problems constitute an important research and ap¬ 
plication domain of tropical mathematics. As an applied mathematical 
discipline that concentrates on the theory and methods of semirings with 
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idempotent addition, tropical (idempotent) mathematics dates back to the 
works of [23, 1 [13, [28] and [25], at least two of which [3, [13] have been 
motivated and illustrated by optimization problems. 

Many subsequent publications that contributed to the development of 
tropical mathematics, including the monographs by and 

a number of contributed papers, were concerned with optimization problems, 
most of which have been drawn from real-world applications in operations 
research and management science. 

Multidimensional tropical optimization problems are generally formu¬ 
lated in the tropical mathematics setting to minimize or maximize linear and 
nonlinear functions defined on vectors over idempotent semifields (semirings 
with multiplicative inverses). The problems may include constraints given 
by linear and nonlinear equalities and inequalities. Many of the problems 
that come from real-world applications and, at the same time, admit so¬ 
lutions in the framework of tropical mathematics have nonlinear objective 
functions defined through multiplicative conjugate transposition of vectors 
(see, e.g., an overview in m. 

There are problems with objective functions that involve the tropical 
algebraic product x~Ax, where A is a given square matrix, x is the 
unknown vector, and x~ is the multiplicative conjugate transpose of x. 
These functions appear in various applications in operations research and 
management science, including problems in project (machine) scheduling 
a El [26] El El [20], location analysis [29] EH [21] , and decision making 
mmm, to name only a few. 

The problem of minimizing the product in question was examined in 
early works [3] IB, 26] by using conventional mathematical techniques. It was 
shown that the minimum in the problem is equal to the tropical spectral 
radius of the matrix A, and attained at the corresponding tropical eigen¬ 
vectors of this matrix. Later, the problem was formulated in the framework 
of tropical mathematics in [3], where a complete solution was proposed by 
reducing to a linear programming problem. Solutions based on tropical 
mathematics were derived in 6, 7J. The results of [71 included an implicit 
description of a complete solution in the form of a vector inequality, and 
provided a computational procedure to solve the inequality. Finally, com¬ 
plete solutions in terms of tropical mathematics to both the problem and its 
generalizations, which have objective functions of an extended form as well 
as additional constraints, were given in pH EH EH]. 

In this paper, we consider a new rather general optimization problem, 
which includes known problems as special cases. We provide a complete solu¬ 
tion in an explicit form on the basis of the approach developed in E3E3ES], 
which introduces an additional variable to represent the values of the objec¬ 
tive function, and then reduces the initial problem to a parametrized vector 
inequality. The minimum of the objective function is evaluated by using the 
solution existence conditions for the inequality. A complete solution to the 
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problem is given by the solutions of the parametrized inequality, provided 
the parameter is set to the minimum value. We discuss the computational 
complexity of the result to show that the solution can be obtained in poly¬ 
nomial time. As a consequence, we propose solutions to new special cases 
of the general problem. 

We apply the results obtained to derive a new complete solution of a 
real-world problem that was drawn from project scheduling (see, e.g., [5] [251 
m for further details on optimal scheduling), and also served to motivate 
the study. We consider a project that consists of activities operating in 
parallel under temporal constraints in various forms, including release dates 
and time windows. For each activity, the flow-time is defined to be the 
time interval between its initiation and completion. The objective is to 
find an optimal schedule that minimizes the maximum flow time over all 
activities. This problem is an extended version of that in mb where a 
less complicated system of temporal constraints is considered. To illustrate 
the solution obtained for the problem, and the computational technique 
implemented by the solution, we present a representative numerical example. 

Note that the problem under examination can be formulated as a linear 
program, and then solved by one of the known solution techniques of linear 
programming. However, these techniques usually take the form of iterative 
algorithms, and do not generally guarantee an explicit closed-form solution. 
Unlike the algorithmic approaches, the proposed solution provides direct re¬ 
sults in a compact vector form suitable for further analysis and practical use. 
Considering, in addition, that the new solution can be calculated in poly¬ 
nomial time, it can certainly serve as a helpful complement and supplement 
to existing solutions. 

The paper is organized as follows. Section [2] includes definitions and no¬ 
tation to be used in the subsequent sections. In Section [3l we present some 
preliminary results, including a binomial identity for matrices and the solu¬ 
tion to linear inequalities. The main result is provided in Section [4l where 
we first offer examples of known optimization problems, then formulate and 
solve a new general problem, discuss the computational complexity of the 
solution, and finally, give solutions to new special cases of the general prob¬ 
lem. Section [5] contains an application of the results in project scheduling, 
and concludes with a numerical example. 

2 Basic definitions, notation and observations 

We start with a short introduction in the context of tropical (idempotent) 
algebra to offer a unified and self-contained framework for the formulation 
and solution of tropical optimization problems in the rest of the paper. 
Below, we follow the notation and results in OH im ng, which form a 
useful basis for the analysis and solution of the problems under study in 
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a systematic manner and in a compact closed form. Further details on 
tropical mathematics at both introductory and advanced levels can be found 

2.1 Idempotent semifield 

An idempotent semifield is an algebraic structure (X, ©, ©, 0,1), where X 
is a nonempty set, © and (g> are binary operations, called addition and 
multiplication, 0 and 1 are distinct elements in X, called zero and one, 
such that (X, ©, 0) is an idempotent commutative monoid, (X \ {©},©,!) 
is an Abelian group, and multiplication distributes over addition. 

The semifield has idempotent addition, which implies that x©x = x for 
each ifX, and invertible multiplication, which allows each nonzero x to 
have its multiplicative inverse x _1 such that x © x -1 = 1. 

Idempotent addition induces a partial order on X such that x < y if 
and only if x © y = y . It follows from this definition that x < x © y and 
y < x © y. Furthermore, both operations © and © are monotone, which 
implies that the inequality x < y yields x ® z < y ® z and x ® z < y ® z 
for all z. The inversion is antitone, which means that the inequality x < y 
results in x ~ 1 > y~ l for nonzero x and y. Finally, the inequality x®y < z 
is equivalent to the two inequalities x < z and y < z. 

It is assumed that the partial order can be extended to a linear one 
to take the semifield as linearly ordered. The relation symbols and the 
optimization objectives are considered below in terms of this order. 

Integer powers are routinely used as shorthand for iterated multiplica¬ 
tion such that x° = 1 and x rn = x © x m ~ 1 for all x S X and integer 
m > 1. Moreover, it is assumed that the equation x rn = a has a solution 
for any a € X and positive integer m, which extends the power notation 
to rational exponents, and thus makes the semifield algebraically complete 
(radicable). In the expressions that follow, the multiplication sign © is 
omitted for brevity. 

Examples of the semifield include R m ax,+ = (RU {—oo}, max, +, —oo, 0) 
and IR m in,x = (R+ U {©oo}, min, x, ©oo, 1), where 1R is the set of real num¬ 
bers and [R + = {x > 0|x G R}, to list only a few. 

The semifield R m ax,+ is equipped with addition and multiplication de¬ 
fined, respectively, as max and +. Furthermore, the number —oo is taken 
as zero, and 0 is as one. Each x G R has the inverse x _1 , which corresponds 
to the opposite number —x in the usual notation. The power x y exists for 
any x,y G R and coincides with the ordinary arithmetic product xy. The 
order defined by idempotent addition is consistent with the conventional 
total order on R. 

In Rmm, x j we have ffi = min, © = x, © = ©oo and 1 = 1. The 
inversion and exponentiation notations have the usual meaning. The relation 
< defines an order that is opposite to the standard linear order on R. 
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2.2 Matrices and vectors 


The set of matrices of m rows and n columns over X is denoted X mxn . A 
matrix with all entries equal to 0 is the zero matrix denoted by 0. A matrix 
is row- (column-) regular, if it has no zero rows (columns). 

Matrix addition and multiplication, and scalar multiplication follow the 
usual rules with the scalar operations © and © in place of the ordinary addi¬ 
tion and multiplication. The above inequalities, which represent properties 
of the scalar operations, are extended entry-wise to matrix inequalities. 

For any matrix A € X mXTl , its transpose is the matrix A 1 £ X nxm . 

The square matrices of order n form the set denoted by X nxn . A square 
matrix having 1 along the diagonal and 0 elsewhere is the identity matrix 
denoted by I. For any square matrix A, the nonnegative integer power is 
defined as A 0 = I and A m = AA m ~ l for all integers m > 1. 

The trace of a matrix A = (a*j) € X nxn is given by 

n 

tr A — CL\\ © * * * © CLnn — © i. 

i =1 

The trace possesses the usual properties given by the equalities 
tr(A © B) = tr A © tr B, tv(AB) = tr(BA), tr(xA) = xtrA, 

which are valid for any matrices A, B £ X nxn and scalar x £ X. 

A matrix with only one column (row) is a column (row) vector. In what 
follows, all vectors are column vectors unless otherwise indicated. The set 
of column vectors of order n is denoted X n . 

A vector is regular if it has only nonzero elements. Let x £ X n be a 
regular vector and A £ X nxn be a row-regular matrix. Then, the result of 
the multiplication Ax is a regular vector. If the matrix A is column-regular, 
then the row vector x T A is regular as well. 

For any nonzero vector x £ X n , its multiplicative conjugate transpose 
is the row vector x~ = (xj~), where xj = x^ 1 if x t ^ 0, and = 0 
otherwise. 

The conjugate transposition exhibits some significant properties to be 
used later. Specifically, if x and y are regular vectors of the same order, 
then the inequality x < y implies x~ > y~ and vice versa. Furthermore, 
for any nonzero vector x , the equality x~x = 1 holds. Finally, if the vector 
x is regular, then the matrix inequality xx~ > I is also valid. 

A scalar A £ X is an eigenvalue of a matrix A £ X nxn , if there exists a 
nonzero vector x £ X n such that Ax = Xx. The maximum eigenvalue is 
referred to as the spectral radius of A, and given by [H [Ml ;25j 

n 

A = tr A © ■ ■ ■ © tr 1//n (A n ) = 0 tr 1 / m (A m ). 

m— 1 
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3 Preliminary results 


We now offer some auxiliary results to be used in the subsequent analy¬ 
sis of optimization problems. We start with binomial identities for square 
matrices, and then describe solutions to linear vector inequalities. 

The inequalities are examined using somewhat different techniques and 
notation by many authors, including |28l |3, [30| GQ HO] - Below, we offer 
solutions given in a compact vector form that provides a natural basis for 
solving the optimization problems in a straightforward and concise manner. 

3.1 Binomial identities 

Let A and B be square matrices of the same order, and m be a positive 
integer. Then, the following binomial identity clearly holds: 

m 

(A© J B) m = 0 0 B io (AB h AB 12 ••• AB ik ) ® B m . 

k =1 -1- ik=m—k 

As an extension of this identity, we derive the following results. First, 
after summation over all m and rearrangement of the output to collect terms 
of like number of cofactors A , we obtain the matrix equality 

mm m 

0(A®S) fc = 0 0 B i °(AB il ■ ■ ■ AB ik ) © 0 B k . (1) 

k =1 k =1 0<io+u4 \-ik<m—k k =1 

Furthermore, by applying the trace and by using its properties, we 
rewrite m in the form of the scalar equality 

mm m 

©tr(.4®B)* = © © ti(AB h ■ ■ ■ AB lk ) ffi 0 tr B k . (2) 

k =1 k =1 0<nH \-ik<m—k k =1 

Both identities m and © are used below to expand matrix expressions 
in evaluating the minimum of the objective function. 

3.2 Linear inequalities 

Suppose that, given a matrix A G x mxn and a regular vector d E X m , the 
problem is to find all vectors a:eX“ that satisfy the inequality 

Ax < d. (3) 

A complete direct solution to the problem under fairly general conditions 
can be found in the following form (see, e.g., Id)- 

Lemma 1. For any column-regular matrix A and regular vector d, all 
solutions to © are given by 

x < (d~A)~. 
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Furthermore, we consider the problem: given a matrix A £ X nxn and a 
vector b £ X n , find all regular vectors a; £ X n that solve the inequality 

Aa? © 6 < x. (4) 

To describe a solution to inequality (JH) in a compact form, we introduce 
a function that maps each matrix A £ X nxn onto the scalar 

Tr( A) = tr A © ■ ■ ■ © tr A n , 

and use the asterate operator (the Kleene star), which takes A to the matrix 

A* = / © A © • • • © A n_1 . 

Presented below is a complete solution proposed in [IS] . 

Theorem 2. For any matrix A and vector b, the following statements hold: 

1. If Tr(A) < 1, then all regular solutions to inequality dH) are given by 
x = A*u, where u is a regidar vector such that u> b. 

2. If Tr(A) > t, then there is no regular solution. 

To conclude this section, we present a solution to a system that combines 
inequality (j4]) with an upper bound on the vector x in the form 

Ax © b < x, 

“ ( 5 ) 

x < a. 

By application of both Lemma |T] and Theorem [2J we arrive at the next 
solution, which is also a direct consequence of the result obtained in [16] for 
a slightly more general system. 

Lemma 3. For any matrix A, vector b and regular vector d, we denote 
A = Tr(A) © d~A*b. Then, the following statements hold: 

1. If A < 1, then all regular solutions to system (JSJ) are given by x = 
A*u, where u is a regular vector such that b < u < (d~A*)~ . 

2. If A > 1, then there is no regular solution. 


4 Solution to optimization problems 

In this section, we consider optimization problems involving the function 
x~Ax, where A is a given matrix, and x is the unknown vector. The un¬ 
constrained minimization of this function is examined by different methods 
in various application contexts HUHU261 [3U6U7] . Complete solutions to some 
constrained problems are proposed in mmm 
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We present examples of both unconstrained and constrained problems, 
and then formulate and solve a new general constrained optimization prob¬ 
lem. As a consequence, we offer solutions for some new special cases of the 
general problem. 

The results are given in the context of an arbitrary idempotent semifield 
in a common form, which can be readily interpreted in terms of particular 
semifields. Specifically, for the semifield K m ax,+, we replace © by max and 
ffi by +, and use the relation symbol < in the usual sense. In the framework 
of IRmin, x j we put © = min and ® = x , and understand the symbol < to 
indicate the order, which is opposite to the standard linear order on [R. 

4.1 Examples of optimization problems 

We start with an unconstrained problem that has the objective function 
written in a basic form. Given a matrix A £ X nxn , consider the problem 
to find regular vectors x £ X n that 

minimize x~Ax, (6) 

A solution to the problem can be provided by several ways (see, e.g., 
DUES DU), and takes the following form. 

Lemma 4. Let A be a matrix with spectral radius A > 0. Then, the 
minimum value in problem ([6]) is equal to A, and all regular solutions are 
given by 

x = (A _1 A)*u, u £ X n . 

Some extensions of problem d6j) were examined in HU E8J El, where 
more general forms of the objective function are considered and/or further 
inequality constraints are added. Specifically, a problem with an extended 
function is solved in m Given a matrix A £ X nxn , vectors p, q £ X n , and 
a scalar r £ X, the problem is to obtain regular x £ X n that 

minimize x~ Ax © x~p ffi q~ x © r. (7) 

A complete direct solution to the problem is as follows. 

Theorem 5. Let A be a matrix with spectral radius A > 0, and q be a 
regular vector. Then, the minimum value in problem © is equal to 


n—1 

fJL = A ffi ® (q- A m p) l ^ m+ ^ ffi r, 

m =0 

and all regular solutions are given by 

x = ( n~ l A)*u , p~ 1 p < u < p(q~(p~ 1 A)*) - . 
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Suppose now that, given matrices A, B G X nxn , and a vector g G X n , 
we need to find regular solutions x G X n to the problem 

minimize x~ Ax , 

( 8 ) 

subject to 5a? © (7 < a?. 

The next complete solution to the problem is provided in [IB]. 

Theorem 6. Let A be a matrix with spectral radius A > 0, and B a matrix 
with Tr (B) < 1. Then, the minimum value in problem © is equal to 

71—1 

/i = A©® 0 tr 1 ^ (AB h ■■■AB Lk ), 

k =1 1<uH- \-ik<n—k 

and all regular solutions are given by 

x = (pT l A © B)*u, u>g. 

Below, we offer a solution to a new problem that combine the objective 
function in (]7|) with the extended set of constraints in Q. 

4.2 New constrained optimization problem 

This section includes a complete solution to a constrained problem, which 
presents an extended version of the problems considered above. We follow 
the approach developed in OS ESDI] to introduce an additional variable, 
which represents the minimum value of the objective function, and then to 
reduce the problem to an inequality, where the new variable plays the role 
of a parameter. 

Suppose that, given matrices A,B £ X nxn , vectors p,q,g,h G X n , and 
a scalar r G X, the problem is to find regular vectors x G X n that 

minimize x~~ Ax © x~ p © q~x © r, 

subject to Bx (B g < x , (9) 

x < h. 

We start with some general remarks and useful notation. It immediately 
follows from Lemma [3] that the inequality constraints in Q have regular 
solutions if and only if the condition Tr(B) ® h~B*g < 1 holds, which is 
itself equivalent to the two conditions Tr(i?) < 1 and h~B*g < 1. 

Clearly, the constraints can be rearranged to provide another represen¬ 
tation of the problem in the form 

minimize x~ Ax © x~p © q~x © r, 

subject to Bx < x, (10) 

g < x < h. 
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To describe the solution in a compact form, we introduce an auxiliary 
notation for large matrix sums. First, we define the matrices Sq = I and 

S k = © AB’ 1 • • • AB ik , k = l,...,n- (11) 

0<iiH- \-ik<n—k 

and note that they satisfy the inequality S k > A k . 

For a different type of sums, we introduce the notation To = B* and 

T k = © B io (AB il ■ ■ ■ AB ik ), k = 1,... ,n — 1. (12) 

0<zo+nH- \-ik<n—k—l 

It is easy to see that the matrices are related by the equality S k+ \ = 
AT k , which is valid for all k = 0,1,... , n — 1. Finally, note that, under the 
condition B = 0, the matrices reduce to S k = A k and T k = A k . 

We are now in a position to offer a complete solution to problem ((9|). 

Theorem 7. Let A be a matrix with spectral radius A, and B be a matrix 
such that Tr (B) < 1. Let p and g be vectors, q and h be regidar vectors, 
and r be a scalar such that h B*g < t and A © ( q~ p) 1 ^ 2 © r > 0. 

Then, the minimum value in problem Q is equal to 

n n —1 

0 = ©tr 1 /*^*) © ©(h'T^) 1 /* 
k=1 k =1 

n —1 n—1 

© Q)(q~T k g © h~T k p) 1 /^ © ©(q-T^) 1 /^) © r , 
k =0 k =0 

and all regular solutions are given by 

x = (1 9~ 1 A®B)*u , 

where u is any regular vector that satisfies the conditions 

6~ l p ®g<u< (( 9~ l q~ © h-){9~ l A © B)*)~. 

Proof. We introduce a parameter to represent the minimum value of the 
objective function, and then reduce the problem to solving a parametrized 
system of linear inequalities. The necessary and sufficient conditions for the 
system to have regular solutions serve to evaluate the parameter, whereas 
the general solution of the system is taken as a complete solution to the 
initial optimization problem. 

Denote by 9 the minimum of the objective function over all regular 
vectors x. Then, all regular solutions to problem ([U]) are determined by the 
system 

x~ Ax © x~p © q x © r < 9, 

Bx © g < x, (13) 

x < h. 
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The first inequality at (fT3l) is equivalent to the four inequalities 

x Ax < 0, x~p < 0, q~x < 0, r <9. (14) 

We use these inequalities to derive a lower bound for 0 and verify that 
9 ^ 0. The first inequality at (1141) and Lemma [4] imply that 9 > x~Ax > A. 
From the next two inequalities and a property of the conjugate transposition, 
we derive 9 2 > q ~xx p > q p, which gives 9 > (q p) 1 / 2 . Since 9 > r as 
well, we finally obtain a lower bound for 9 in the form 

6 > A © (q'p ) 1 / 2 © r, (15) 

where the right-hand side is nonzero by the conditions of the theorem. 

We can now multiply the first two inequalities at (1141) by 0 _1 , and then 
apply Lemma [3] to the first three. As a result, we have the inequalities 

9~ l Ax < x, 9~ 1 p < x, x < 9q. 

As the next step, we combine these inequalities with those in the system 
at m Specifically, the first two inequalities together with Bx © g < x 
give the inequality (9~ l A © B)x © 9~ 1 p © g < x. 

In addition, we take the inequalities x <9q and x < h, and put them 
into the forms x~ > 9~ 1 q~ and x~ > h ~. The last two inequalities are 
combined into one, which is then rewritten to give x < ( 9~ 1 q~ © h~)~. 

By coupling the obtained inequalities, we represent system m as 

(0 _1 A © B)x © 0 -1 p © g < x, 

, 16 
x < (9 l q © h ) . 

Considering that system (1161) has the form of we can apply Lemma[3] 
to examine this system. By the lemma, the necessary and sufficient condition 
for (1161) to have regular solutions takes the form 

Trie- 1 A © B) © (0-V © h~){9~ l A © B)*{9~ 1 p © g) < 1 . 

To solve this inequality with respect to the parameter 0, we put it in a 
more convenient form by expanding the left-hand side in powers of 0 . 

As a starting point, we examine the matrix asterate 

71—1 71—1 

[9~ l A © B)* = ®(0 _1 A © B) k = I © ®(0 _1 A © B) k . 

k =0 k =1 

After application of © to the second term, we rearrange the expression 
to collect terms with the same power of 0 , and then use m to write 


{e~ 1 A®B)* 


n —1 


© © 

k=l 0<7q+HH k —1 


n— 1 

9~ k B l0 {AB n ■ • • AB ik ) © ® B k 

k =0 


71—1 71—1 

= ®0- fc T,©T o = ®0- fc T,. 

k =1 fc=0 
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By using ©, (ED, G2D and properties of the trace function, we also 
have 

n 

Tr(0 _1 A ® B) = 0 tr^ 1 A © B) k 

k =1 

n n 

= © © 9~ k tv(AB h ■■■ AB ik )®Q)tr(B k ) 

k =1 0<uH- \-ik<n—k k =1 

n 

= ®0- fe tr(5 fc )0Tr( J B). 
fc=i 

Substitution of these results into the condition for regular solutions yields 

n n —1 

© 0~ k tr (S k ) © ® e-^e-'q- © h-)T k (6- l p © g) © Tr(B) < 1. 
k =1 fc=0 

Since Tr(_B) < 1 by the conditions of the theorem, the term Tr(£>) 
does not affect the solution of the inequality, and hence can be omitted. 
The remaining inequality is equivalent to the system of inequalities 

9~ k tr (Sk) < 1, k = 1,..., n; 

9~ k {0~ l q~ © h~)T k (0~ l p © g) < t, k = 0,1,..., n — 1; 

which can be further split into the system 

9~ k tr (Sk) < 1, k = 1 ,..., n; 

9~ k h-T k g < 1, 

9~ k ~ 1 (q T k g © h~T k p) < 1, 

9~ k ~ 2 q~T k p < 1, k = 0,1,..., n — 1. 

Note that h T 0 g = h~B*g < 1 by the conditions of the theorem, and 
thus the second inequality in the system is valid at k = 0 for all 9 > 0. 

By solving the inequalities, we have 

9>tr 1/k (S k ), fc = l,...,n; 

9>(h-T k g) 1 / k , k = l,...,n —1; 

9> {q-T k g®h-T k p) l ^ k+l \ 

9>{q-T k p) 1 ^ k+2 \ fc = 0,l,...,n-l. 

The obtained solutions can be combined into one equivalent inequality 

n n —1 

9 > ©tr ^(Sk) © ©(/i-Tfcff) 1 ^ 

k =1 /c=l 

n—1 n—1 

© ©(q-T^ © © ® (qrTfcp) 1 /^ 2 ). 

fc=0 /c=0 
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We have to couple the lower bound given by (11511 with that defined by 
the last inequality. It is not difficult to verify that the right-hand side of this 
inequality already takes account of the terms A and (q^p) 1 ^ 2 presented in 
(USD- Indeed, considering that S k > A k , we have 

n n 

©tr 1 / fc (S fc )>©tr 1 / fc (A fc ) = A. 

k= 1 k= 1 

Moreover, since To = B* > I, it is easy to see that 

71—1 

®(q-T fc p) 1 /(fc+ 2 ) > (q-Top ) 1 / 2 > (<T;p) 1/2 . 

k=0 

By combining all lower bounds obtained for 9 , we arrive at the inequality 

71 71—1 

0 > ® trV^Sfc) © ®(^-T fc q) 1 / fc 
fc=l fc=l 

71—1 71—1 

0 ®(<TT fc q © fe-Tfcp) 1 /^ 1 ) © ® {q-T k pf/^ © r. 

/c=0 k =0 

Since 9 is assumed to be the minimal value of the objective function, 
this inequality must hold as an equality, which yields the desired minimum. 

Finally, we take the minimum value of 9, and then apply Lemma [3] to 
obtain all solutions of the system at (1161) in the form 

x = ( 9~ l A © B)*u, 9~ l p © g <u < {(9~ l q~ © h~)(9~ 1 A © B)*)~. 

Because the solution obtained is also a complete solution of the initial 
optimization problem, this ends the proof of the theorem. □ □ 

We conclude this section with a brief discussion of the computational 
complexity of the solution obtained to see that it is polynomial in the di¬ 
mension n. Indeed, this complexity is determined by the complexity of 
computing the minimum value 9 , as the other components of the solution 
are given by a finite number of matrix and vector operations, and thus ob¬ 
viously take no more than polynomial time. 

Furthermore, it directly follows from the expression for 9 that, if the 
evaluation of the matrix sequences S\ : .... S n and To,..., T n -\ has poly¬ 
nomial complexity, then so has that of 6. Considering that S^+i = AT for 
all k = 0 ,..., n — 1, we need to verify that T& can be obtained in polynomial 
time. 

To describe a polynomial scheme of calculating TJt, we first write 

n—k—1 

T k = ® Qku Qki = © 

l = 1 -b 'i'k = L 
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where Qki is the sum of all matrix products that are comprised of k factors 
equal to A and l factors equal to B , with Q^o = A k , Q 0 i = B l and 
Qoo = I In this case, the evaluation of the matrices To,..., T n _i reduces to 
computing the matrices Qki for all k = 0 ,..., n — 1 and l = 0 ,..., n — k — 1 . 

Furthermore, we note that the recurrent relation Qkl — -AQk—1,1 © 
holds for all k,l = 1 , 2 ,... It is clear that this relation offers a 
natural way to obtain successively all matrices Qki , using two matrix mul¬ 
tiplications and one matrix addition per matrix. Since the overall number 
of matrices involved in computation is 1 + 2 + ■ ■ ■ + (n — 1 ) = n{n — l)/ 2 , 
the computation of all matrices T^ requires polynomial time, and thus the 
entire solution has polynomial complexity. 

4.3 Some special cases 

As direct consequences of the result obtained, we now find solutions to 
special cases of problems Q and cnD with reduced sets of constraints. To 
begin with, eliminate the first constraint in (| 10 l) and consider the problem 

minimize x~ Ax © x~p © q~x © r, 

( 17 ) 

subject to g < x < h. 

Clearly, the solution to this problem can be derived from that of (11011 
by setting B = 0. Under this condition, we have Sk = A k and T& = A k , 
whereas the solution is described as follows. 

Corollary 8 . Let A be a matrix with spectral radius A. Let p and g be 
vectors, q and h be regular vectors, and r be a scalar such that h~g < 1 
and A © (q~p) 1 / 2 © r > 0. Then, the minimum value in problem m is 
equal to 


n —1 n —1 

6 = A © Q)(h~A k g) l / k © ® {q~ A k g © hT A fc p) V(fc+1) 

k =1 k =0 


77,-1 

© ® (q~ A k p) l/( - k+2) © r, 

k =0 


and all regular solutions are given by 


x = (6 1 A)*u , 9 l p © g < u < ((9 1 q © h )(9 1 A)*) . 


Furthermore, we consider another special case of (110(1 . which takes the 
form 

minimize x~ Ax © x~p © q~x © r, 

( 18 ) 

subject to Bx < x. 

After slight modification of the proof of Theorem [71 we arrive at the next 
result, which can also be obtained directly by putting g = 0 and h~ = 0 T 
in the solution of problem ( 1101 ) . 
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Corollary 9. Let A be a matrix with spectral radius X, and B be a matrix 
such that Tr (B) < 1. Let p be a vector, q be a regular vector, and r be a 
scalar such that X®{q~p) 1 ^ 2 ®r > 0. Then, the minimum value in problem 
m is equal to 


n n— 1 

0 = 0 tr V k (S k ) © 0 (q-Tkp) 1 /^ © r , 
k=1 k=0 

and all regular solutions are given by 

x = {0~ l A © B)*u , e^p <u< e{q~(e~ 1 A 0 B)*)~. 

Finally, note that eliminating both inequality constraints in ([U]) leads to 
the same solution as that provided by Theorem [5j 

5 Application to project scheduling 

We now apply the result obtained to solve an example problem, which is 
drawn from project scheduling EE3U27] and serves to motivate and illus¬ 
trate the study. 

Consider a project consisting of a set of activities that are performed 
in parallel under various temporal constraints given by precedence relation¬ 
ships, release times and time windows. The precedence relationships are 
defined for each pair of activities and include the start-finish constraints on 
the minimum allowed time lag between the initiation of one activity and 
completion of another, and the start-start constraints on the minimum lag 
between the initiations of the activities. Once an activity starts, it con¬ 
tinues to its completion, and no interruption is allowed. The activities are 
completed as soon as possible under the start-finish constraints. 

The release time constraints take the form of release dates and release 
deadlines to specify that the activities cannot be initiated, respectively, be¬ 
fore and after prescribed times. The time windows are given by lower and 
upper boundaries, and determine the minimum time slots preallocated to 
each activity. The activities have to occupy their time windows entirely. If 
the initiation time of an activity falls to the right of the lower boundary of 
its window, this time is adjusted by shifting to this boundary. In a similar 
way, the completion time is set to the upper boundary if it appears to the 
left of this boundary. 

Each activity in the project has its flow-time defined as the duration of 
the interval between the adjusted initiation and completion times. A sched¬ 
ule is optimal if it minimizes the maximum flow-times over all activities. 
The problem of interest is to find the initiation and completion times of the 
activities to provide an optimal schedule subject to the temporal constraints 
described above. 
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5.1 Representation and solution of scheduling problem 

Suppose a project involves n activities. For each activity i = 1 ,,n, let 
Xi be the initiation and y* the completion time. We denote the minimum 
possible time lags between the initiation of activity j = 1 and the 

completion of i by a t j , and between the initiations of j and i by b tJ . If a 
time lag is not specified for a pair of activities, we set it to —oo. 

The start-finish constraints yield the equalities 

Vi = max(uji + xi,..., a in + x n ), i = 1,..., ra; 
whereas the start-start constraints lead to the inequalities 

Xi > max(6ji + x\ ,b in + x n ), i = l,...,n. 

Let gt and h % be, respectively, the possible earliest and latest initiation 
times. The release date and release deadline constraints are given by the 
inequalities 

9i < Xi < hi, i = l,...,n. 

Then, we denote the lower and upper boundaries of the minimum time 
window for activity i by qi and pt, respectively. Let be the adjusted 
initiation time and ti the adjusted completion time of the activity. Since 
the time window must be fully occupied, we have 

Si = mm(xi,qi) =-m&x(-Xi,-qi), U = max(j/ i ,}) i ) ) i = 

Finally, the maximum flow-time over all activities is given by 
max(b - si,... ,t n - s n ). 

We are now in a position to represent the optimal scheduling problem of 
interest as that of finding Xi, yi, s* and ti for all i = 1, ... ,n to 

minimize max (tj — s,), 

l<i<n 

subject to Si = — max(—Xj, -%), U = ma x(yi,pi), 

Hi = max (aij +Xj), Xi > max (bij + Xj), 

9i <Xi<hi, i = 1,. ., ,n. 

It is not difficult to see that this problem can be represented and solved 
within the framework of linear programming, which generally offers algo¬ 
rithmic solutions rather than a direct complete solution in an explicit form. 

To obtain a direct solution, we place the problem in the context of trop¬ 
ical mathematics. Considering that the problem is formulated only in terms 
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of the operations of maximum, ordinary addition, and additive inversion, 
we can rewrite it in the setting of the semifield K max ,+ as follows: 


minimize 
subject to 


n 


(D s i 


i= 1 


Si = (xr 1 eg" 1 )' 1 , 

ti = Vi® Pi: 

n 

n 

© 

ii 

s 

© 

Al 

3 = 1 

3 =1 

9i<Xi< hi, 

i = 1 ,...,n 


Furthermore, we put the problem into a compact vector form. We in¬ 
troduce the matrix-vector notation 


A = ( cHj ), B = ( bij ), x = (xi), y = (y*), g = (, gi ), h = (hi), 

and write the start-finish, start-start and release time constraints as 

y = Ax, x > Bx, g < x < h. 

To take into account the time window boundaries and adjusted times, 
we use the vector notation 


s = (si), t = ( ti ), p = ( Pi ), q = ( qi )• 

The vectors of adjusted initiation and completion times take the form 
s = (x~ © q~)~, t = y®p. 

The optimal scheduling problem to minimize the maximum flow-time 
subject to the temporal constraints under consideration now becomes 

minimize s~t, 

subject to s - = x~ © q~, t = y®p, 

Ax = y, Bx < x, 

g < x < h. 

Note that, in the context of scheduling problems, it is natural to consider 
the matrix A as column-regular matrix, and the vectors p, q and h as 
regular. 

A complete solution to the problem is given by the next result. 

Theorem 10. Let A be a column-regular matrix, and B be a matrix such 
that Tr(£?) < 1. Let p, q and h be regular vectors and g be a vector such 
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that h B*g < 1. Then, the minimum flow-time in problem CED is equal to 

n n—1 n 

0 = 0 t r^(Sfc) © ©(/TT^) 1 /* © Q(q~S k g) l ' k 

k =1 k =1 k =1 

n—1 n 

© ©(fr-Tfcp) 1 /^ 1 ) © ©(g-^fcp) 1 /^ 1 ), (20) 

k =0 fc=0 

and, the vectors of initiation and completion times are given by 

x = ( 8 ~ 1 A@B)*u, y = A{ 6 ~ 1 A © B)*u, (21) 

8 = ((( 6 ~ 1 A®B)*u)- ®q~y, t = A{e~ 1 A®BYu®p, ( 22 ) 

where u is any vector that satisfies the conditions 

e~ l p ®g<u< (0r V A © h~)(9~ 1 A © B)*)~ . (23) 

Proof. First, we eliminate the vectors s and t from problem (HID by repre¬ 
senting the objective function as 

s~t = ( x~ © q~)(y © p) = x~y © q~y © x~p © q p. 

Furthermore, we substitute y = Ax to reduce (fl9l) to the problem 

minimize x~ Ax © q~ Ax © x~p © q p, 
subject to Bx < x, 

g < x < h, 

which has the form of (1101) . where q~ is replaced by q A and r by q~p. 

To apply Theorem 0 we note that, under the given conditions, the con¬ 
ditions of the theorem are satisfied as well. Specifically, since both vectors 
p and q are regular, we have r = q~p > 0, and thus provide the last 
condition of Theorem [71 

Next, we refine the expression for 8 by applying the identity AT k = 
S k+ \ , which is valid for all k = 0 , ..., n — 1. 

After some rearrangement of sums, we arrive at (1201) . Both the repre¬ 
sentation for x at (1211) and the condition on u at (1231) are directly obtained 
from Theorem [7J The other expressions in m and (1221) are immediate 
consequences. □ □ 

As before, the solutions to special cases without constraints are readily 
derived from the general solution offered by Theorem [TUI Specifically, we 
eliminate the boundary constraint g < x < h by setting g = 0 and h~ = 
0 7 , and/or the linear inequality constraint with matrix in the form Bx < x 
by setting B = 0, which further yields the substitutions S k = A k and 
T k = A k . 
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5.2 Numerical example 

To provide a clear illustration of the above result and of the computational 
technique, we solve in detail a simple low-dimensional problem. Even though 
the example under consideration is somewhat artificial, it well demonstrates 
the applicability of the solution to real-world problems of higher dimension. 

Let us examine a project that involves n = 3 activities under constraints 
given by the matrices 



We start with the verification of the existence conditions for regular 
solutions in Theorem 1 101 First note that the matrix A is obviously column- 
regular. In what follows, we need the powers of the matrix A, which have 
the form 

/ 8 4 1 \ / 12 8 5 \ 

A 2 = 6 6 4 , A 3 = 10 9 7 . 

\ 5 4 6 / \ 9 7 9 / 

Then, we take the matrix B and calculate 


/ 0 

—oo 

1 \ 

/ 0 

-l 

1 \ 



-1 

1 ’ 

B 3 = [ 0 

0 

2 ’ 

Tr(£>) = 0 

\ —oo 

-2 

0 / 

V -i 

—oo 

0 / 



Furthermore, we successively obtain 
( 0 -1 1 \ 

B* = 1 0 2, h~B* = ( -2 -3 -1 ) , h~B*g = 0. 

V -1 -2 0 ) 


Since Tr (B) = h B*g = 0 , where 0 = 1 , we conclude that the con¬ 
ditions of Theorem m are fulfilled, and thus the problem under study has 
regular solutions. 

As the next step, we find the minimum value 0 by application of ( 1201 ) . 
The evaluation of 9 involves the matrices 


S 0 = I. = A © AB © AB 2 , S 2 = A 2 © ABA® A 2 B. S 3 = A 3 , 

T 0 = B *, Tj = A © AB © BA, T 2 = A 2 . 
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To obtain S \, S -2 and T \, we calculate the matrices 


/ 0 3 5 \ / 2 2 4 \ 

AS = 315 , BA= 4 3 5 , 

\ 2 0 3 / \ 3 -1 -oo / 

and then the matrices 

/ 4 —1 5 \ / 6 6 8 \ / 4 7 9 \ 

AS 2 = 4 2 4, ABA = 7 6 8 , A 2 # = 6 5 8 . 

\ 2 13/ \ 6 4 6 / V 546 / 

After substitution of these matrices, we have 

/ 4 3 5 \ / 8 7 9 \ / 4 3 5 \ 

Si= 4 3 5 , £> 2 =j 7 6 8 1 , Ti= 435 . 

\ 2 1 3 / \ 6 4 6 / \ 3 1 3 / 

Based on the results obtained, we calculate the sum 

3 

®tr 1 / fe (S fe ) = 4. 

k =1 

To evaluate the remaining sums, we first find the vectors 
hTTo = ( -2 -3 -1 ) , B^Ti = ( 2 1 3 ) , h To = ( 6 3 3 ) , 
and then obtain 

hT\g = 4, hTT 2 g = 6 , h~Top = 4, h~T\p = 8 , h,~T 2 p = 10. 

With these results, we get another two sums 

(g^-Tfcs) 1 /* = (g^-Tfcp) 1 /^ 1 ) = 4. 

fc=l fc=0 

Furthermore, we obtain the vectors 

<TSo=( -3 -2 -1 ), <T«i = ( 2 13). 

q~S 2 = ( 5 4 6 ) , = ( 9 7 8 ) , 

and then calculate 

qS L g = 4, qS 2 g = 7, = 9, 

q~SoP = 4, q~ S\p = 8 , g _ S 2 p = 11, g~S 3 p = 13. 
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Finally, we use the above results to find the last two sums 


®(<rW A = ®(<r.W /(fe+1) = 4. 

k =1 fc=0 

By combining all sums according to (1201) . we have 

6 = 4. 

To describe the solution set defined by (12T1) and (1231) . we first obtain 
e~ 1 q~A=(- 3 -3 -2), A © h~ = ( -2 -3 -2). 

We calculate the matrices 


0~ 1 A®B = 


and then find 


0 -1 1 
0-1 2 
-1 -3 -1 


{ 0 ~ 1 A®B)* = 


0 -1 1 

(0~ l A © B) 2 = | 1 -1 1 

- 1-2 0 

0 -1 1 
1 0 2 
-1 -2 0 


With (12T1) . all solutions x = ( x\,X 2 ,xs) T to the problem are given by 
x = {6~ X A © B)*u, u\ < u < U 2 , 
where the bounds for the vector u = (ui,tt 2 ,rt 3 ) T in (f251) are defined as 


U\ = 9 p®g = 0 , u 2 = ((9 q A®h ){9 L A®B)*) = 3 


1 


1 


Note that the columns in the matrix (9 l A © B)* are equal up to con¬ 
stant factors, and therefore, this matrix can be represented as 


0 -1 1 

1 0 2 | = 

-1 -2 0 


1 

2 

0 


( -1 -2 0 ). 


We introduce a new scalar variable 

v = ( -1 -2 0 ) u, 

and rewrite the solution in the form 


x = \ 2 | v, 


Vl < V < V2, 
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where the lower and upper bounds on v are given by 


v\ = ( —1 —2 0 ) U\ = 1, V ‘2 = ( — 1 —2 0 ) U 2 = 1. 

Since both bounds coincide, we have the single vector of initiation time 



Finally, using formulas (1211) and (1221) gives the vector of completion time 
and the vectors of adjusted initiation and completion times 
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